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This paper deals with the optimal streaky perturbations (which maximize the per- 
turbed energy growth) in a wedge flow boundary layer. These three dimensional 
perturbations are governed by a system of linearized boundary layer equations 
around the Falkner-Skan base flow. Based on an asymptotic analysis of this sys- 
tem near the free stream and the leading edge singularity, we show that for acute 
wedge semi-angle, all solutions converge after a streamwise transient to a single 
streamwise-growing solution of the linearized equations, whose initial condition 
near the leading edge is given by an eigenvalue problem first formulated in this 
context by Tumin (2001). Such a solution may be regarded as a streamwise evolv- 
ing most unstable streaky mode, in analogy with the usual eigenmodes in strictly 
parallel flows, and shows an approximate self-similarity, which was partially known 
and is completed in this paper. An important consequence of this result is that the 
optimization procedure based on the adjoint equations heretofore used to define 
optimal streaks is not necessary. Instead, a simple low-dimensional optimization 
process is proposed and used to obtain optimal streaks. Comparison with previous 
results by Tumin and Ashpis (2003) shows an excellent agreement. The unstable 
streaky mode exhibits transient growth if the wedge semi-angle is smaller than 
a critical value that is slightly larger than tt/6, and decays otherwise. Thus the 
cases of right and obtuse wedge semi-angles exhibit less practical interest, but they 
show a qualitatively different behavior, which is briefly described to complete the 
analysis. 
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I. INTRODUCTION 

Internal streaks (or Klebanoff modes^) of a two-dimensional laminar boundary layer 
denote three-dimensional disturbances that evolve slowly streamwise and show a fast os- 
cillation (with a wavelength comparable to the boundary layer thickness); see Refs.-~- for 
some recent experiments. The linear stability analysis^ demonstrated that three dimen- 
sional perturbations with streaky structure are supported in inviscid shear flows, whose 
energy grows algebraically in time. This initial growth together with a subsequent decay 
due to viscous dissipation is known as transient growth-. Streaks can be forced either 
internally, from obstacles near the leading edge or externally, from perturbations in the 
free stream, and generally interact with the transversal Tollmien-Schlichting modes, either 
enhancing^ or delaying^ transition to turbulence, depending on the streak amplitude. In 
the former case, the effect is known as bypass transition^. 

Luchiniii proposed in the scope of the Blasius boundary layer an analytical description 
of streaky perturbations with a small (compared to 1/5, where S is the boundary layer 
thickness) spanwise wavenumber, which is equivalent to the limit of small distance to the 
leading edge. The extension of this analysis to the Falkner-Skan^ profile was made by 
Turning, who derived an eigenvalue problem describing the velocity components of the 
streaky perturbation. He computed the largest eigenvalue of this problem in terms of the 
wedge angle and found that the unbounded growth is suppressed as the angle exceeds a 
threshold value, which is slightly larger than tt/3. 

An important issue related to streaky perturbations is to find the initial disturbance 
that maximizes the perturbed energy streamwise growth; the resulting streamwise evolv- 
ing streak is known as the optimal streak. For strictly parallel flows, an explicit computa- 
tion of disturbances resulting in the maximum transient growth is possible by optimizing 
over the eigenmodes of the Orr-Sommerfeld operator, as has been done for a number of 
particular parallel flows, both in the contexts of temporal^ - — and spatial stability^. 

For the non-parallel case, in the scope of Blasius boundary layer, Luchini^ and Ander- 
sson et al.— calculated optimal streaks using a method based on the adjoint formulation, 
which is commonly employed in optimal-control problems for distributed parameter sys- 
tems. Using this method, Tumin and Ashpis^ computed the optimal streaks in a Falkner- 
Skan boundary layer and showed the effects of the spanwise wave number and the wedge 
angle on transient growth. In particular, they found that an adverse pressure gradient 
increases the amplification whereas a favorable presure gradient has the opposite effect 
(in accordance with the previous asymptotic result by Tumin^). Levin and Henningson^ 
also studied the optimal disturbance in a Falkner-Skan base flow, obtaining results that 
were consistent with those by Tumin and Ashpis^. 

In a recent paper—, two of the authors analyzed streaks in a Blasius boundary layer 
and showed that, after an initial transient, they approach a unique (up to a constant 
factor) 'mode', which was called the unstable streaky mode. This mode is calculated 
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from the streamwise evolving linearized equations with well-defined initial conditions near 
the leading edge, which are given by the first eigenmode of the eigenvalue problem first 
formulated by Luchint^k The unstable streaky mode provides the optimal streaks when 
the initial conditions are taken sufficiently close to the leading edge (say, x < 10~ 3 ), 
making it unnecessary the optimization process in this limit. In fact, the asymptotic value 
of the optimal spanwise wavenumber was found to be 0.484, which was slightly different to 
its counterpart (0.45) calculated by Luchini and Andersson et al. The latter corresponds 
to an only moderately small value of the initial streamwise coordinate, x ~ 0.01. The 
analysis in Ref.— was based on three main ingredients: 

1. The exact self-similarity of the boundary layer, which allows to eliminate the 
wavenumber from the formulation. 

2. The modal structure near the leading edge, which allows to expand any initial 
condition into a complete system of eigenmodes. 

3. An asymptotic analysis of the behavior of the solutions near the free stream. This 
was essential to both understanding the structure of the solutions and formulating 
a quite efficient numerical scheme, which provided streaks with initial conditions at 
a section extremely close to the leading edge. 

In addition, the approximate self-similarity already detected experimentally in Ref.- and 
confirmed in Ref.— (namely, the wall-normal profile of the streamwise velocity compo- 
nent of optimal streaks, rescaled with its maximum, remains constant streamwise) was 
completed noting that the wall-normal profile of a certain combination of the cross flow 
velocity components shows the same property as the streamwise velocity component. A 
part of the present paper is an extension of this previous work to the Falkner-Skan bound- 
ary layer. In addition, advantage will be taken of the low dimensional nature of streaks 
to develop a quite efficient and simple method to calculate optimal streaks. This method 
consists in optimizing the perturbed energy gain in a low-dimensional solution manifold 
spanned by a few solutions of the streamwise evolving problem. 

The remaining of the paper starts with the formulation of the problem, in section HT| 
where the asymptotic behaviors near the free stream and the leading edge are analyzed, in 
subsections III Al and III Bl respectively. The streamwise evolution of streaks is considered in 
section liTTl where a quite efficient numerical scheme is presented that relies on the previous 
asymptotic results; the modal nature of streaks is considered in subsection IIII Al Optimal 
streaks are studied in section IIVI where the above mentioned optimization method is 
presented, in subsection IIV AL The body of the paper deals with acute wedge semi- 
angles; the cases of right and obtuse semi-angles are considered in an appendix, at the 
end of the paper. The paper ends with some concluding remarks, in section [V] 
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FIG. 1. The symmetric flow past a wedge with angle ir(3, in the cases (a) /3 = (Blasius flow), 
(b) < p < 1, (c) /3 = 1 (stagnation flow), and (d) 1 < /3 < 2. 

II. GENERAL FORMULATION 

We consider the high-Reynolds-number flow incidenting in a wedge with semi-angle 
Ptt/2 (FigJTJ), where the parameter /3 is named after Hartree^. In particular, we study the 
linear stability of the resulting boundary layer flow, under streaky perturbations that ex- 
hibit a spanwise period comparable to the thickness of the boundary layer, 5* = L*/\/Re, 
where L* is the wall wise length of the portion of the wedge under consideration and 
Re = u*L*/u ^> 1 is the Reynolds number based on L* and the free stream velocity 
u* . Nondimensionalization is made according to the usual boundary layer approxima- 
tion, using the following units: L* and u* for the streamise spatial coordinate x and 
velocity u, respectively, and 5* and u* / \/Re for the wall-normal and spanwise coordinates 
(y,z) and velocities (v,w), respectively; the pressure is scaled with p* (u*) 2 / y/Re. The 
two-dimensional, symmetric base flow, (u, v,w,p) = Vf,, 0,Pb), is the Falkner-Skan™ 
similarity solution 

/'■(c) + (p - i)cno 



u b (x,C)=x*-f>F'(t), v b {x,C) = - 
where ( is the selfsimilar, wall normal coordinate 



(2 - 0)g{x) 



1-/3 



C = y/g(x), with g(x) = x 2 -e, 
and the streamfunction F is the solution of the Falkner-Skan equation 
F" + \FF" + (3(1 — F' 2 )]/ (2 — /3) = in < C < oo, 
F(0) = F'(0) = 0, F'(oo) = 1. 



(1) 
(2) 

(3) 
(4) 
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The Falkner-Skan solution includes the cases (FigJT]) of flat plate boundary layer (3 = 0, 
flow past an acute wedge (0 < (3 < 1), stagnation flow towards a flat plate ((3 = 1), and 
flow in an obtuse wedge (1 < /3 < 2); negative values of (3 involve unphysical reverse 
flow and are only of academic interest. Note that the scaled boundary layer thickness is 
proportional to g(x), which invoking ([2} means that it grows streamwise if /3 < 1 and 
decays if /? > 1, while at (3 = 1 the boundary layer flow is parallel. The remaining of the 
paper will concentrate in the case of acute wedge angle. The case (3 > 1 does not promote 
transient growth and will be briefly analyzed in the Appendix, at the end of the paper. 

Perturbations are considered linearizing around the base flow and decomposing in 
normal modes as 

(u - u b ,v - v b ,w,p-p b )= [U x,C , , , — 7= e , (5) 

V x x xVRe / 

to obtain the following linearized boundary layer (LBL) equations 

xd x U = ~_p 9 < U - d < V ~ a 9{x)W, (6) 

xf%u = da u + JL^ - gg + eg - a^+g - _ F% (7) 

((3 + 1)F' + (g - ljCg - (2 - /3)a 2 ^(x) 2 T/ 
+ 2^ V - 9 ? P, (8) 

xF'd x W = d a W + j^d ( W + 1 2 _ P | 91 J ^ - ^(x)P, (9) 

£7 = V = W = at C = and oo, P = at C = °o, (10) 

where £ is the self similar wall normal coordinate defined in eq.([2]) and d x , d y , ■ ■ ■ denote 
hereafter partial derivatives. Note that the boundary conditions are no-slip at the plate 
and vanishing at the infinity in the wall-normal direction. 

If j3 = 1, then g(x) = 1 and the right hand sides of eqs.(EH9]) are independent of x. 
Otherwise, the spanwise wavenumber can be eliminated from the formulation using the 
scaling 

(x,y,u b ,v b ) = (a^,ay,a^u b ,v b /a), (U,V,W,P) = a^(U,V,W,P). (11) 

In addition, as in Ref.— , we use the new cross flow variable 

H = V + W, (12) 

which anticipating results in subsection III Al below, converges quite fast to zero as ( — > 
oo. Thus, we substitute the spanwise velocity component by this new variable in the 
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formulation. Using all these, eqs.(1B|)- (ITU|) are rewritten as 

idJJ = JZjiW - 9(V - 9(i)V + g(x)H, (13) 

+ W + DP + W - i)(F" - (2 - P)g(i) 2 y _ 8(Pi (15) 

+ + - DCF- - + r - (2 - ^ . _ ^ - _ ^ . (16) 

U = V = H = at C = and oo, P = at ( = °o. (17) 

These equations will referred to below as the modified linear boundary layer (MLBL) 
equations. 



A. Asymptotic behavior near the free stream (( >> 1) for < 1 

As £ — > oo, the streamfunction of the Falkner-Skan base flow behaves as 

F(()~(-ap + 0(e-«- a ^/ 2 ), (18) 

where the constant ap depends on Hartree parameter 0. Thus both F' — 1 and F" decay 
to zero exponentially fast and the LBL equations ((6])-(l9)) can be greatly simplified. In 
particular, the streamwise momentum equation (j7|) becomes 

x^f/ = <9 CC C/ + Yzrfd(U — U. (19) 

This equation is unforced, which means that U — 0. Using this, the remaining LBL 
equations dHJ)-© simplify to 

d c V = ag(x)W, (20) 
*W = + + ^f^ V - *P, (2D 

xd x W = d cc W + Y^-d c W + l-(2-y(a;)V - ^(x)P. (22) 
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These equations can be solved in closed form as follows. We first eliminate the pressure 
manipulating (I2ip and ff22]) as usually, to obtain the following equation for the streamwise 
vorticity Q = d^W — ag(x)V 

a^fi = <9 cc fi + - _ £ d c il H — — Q. (23) 

This equation is also unforced, which means that f2 = 0, namely 

d(W = ag(x)V. (24) 
Excluding divergent behaviors as ( — > oo, eqs. fl2"Ul) and fl2"4l) yield 

V = -W = V 00 e~ a9{x){c - a ?\ (25) 

where is a function of x that remains undetermined but will not be necessary below. 
The pressure P is readily obtained (and seen to behave as (^e~ a9 ^ x ' , ^~ a ^) substituting (1251) 
into ([22]). 

1-/3 

Since (3 < 1, g(x) = x 2 -p is small at small x and the convergence of V, W, and P to 
the final free stream state U = V = W = P = 0is quite slow (see Figj3] below), which 
explains the difficulties encountered in former numerical treatments 2 ^— of (l6l)- (fT0l) that 
did not take into account this behavior. The streamwise velocity component U instead, 
converges quite fast to zero (see Fig J3] below) , as e~^~ a ^ 2 / 2 . Equation ( l25|) shows that 
the same happens with the quantity (cf ( fT2l ) 



H = V + W, (26) 

as anticipated above. Note using eqs.(l20l). ( 1241) . and ( |25l) that H coincides with the 
streamwise vorticity at large £. 



B. Behavior near leading edge (x 1) in the case f3 < 1 

Assuming that the solution behaves as a power of x at small x, the relevant behavior is 
given by an eigenvalue problem first formulated and solved by Tumin^ 2 -, whose formulation 
is re-interpreted here. This will be done taking into account the asymptotic behavior as 
C — > oo encountered in last sub-section, which in conjunction with the Tumin scaling 
suggests the ansatz 

(U, V, H, P) ~ x l - x {U, V, H/g(x), p) e -'™- a P\ (27) 

1-/3 

Substituting this into the MLBL-equations and neglecting 0{g(x)) = 0(x 2 -p )-terms 
yields the following eigenvalue problem 

(1-\)U+ %—i^ U / + V'-H = 0, (28) 
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FIG. 2. Left: The two smallest Tumin eigenvalues (subscript T) and the smallest Chen-Libby 
eigenvalue (subscript CL) in terms of ft, in the range < ft < 1; the first Tumin eigenvalue is 
smaller than one if < ft < ft c = 0.362. Plain circles for ft = 1 come from calculations in the 
Appendix. Right: U, V, H (solid lines), and Uq = — \H\ max U /\U\ m ax. (dashed lines) for the first 
Tumin mode and ft = 0.1,0.2,0.353 (that value considered in Ref.—), and 0.5; arrows indicate 
increasing ft. 



u" 



F & (1-^+1(2-^-1]^ _ F „ - = 



2- ft 



2- ft 



y„ + J_y, + (ft-l)[F+(2ft-l)(F>+(ft-l)C 2 F"] u 



2- ft 



+ 



(2 -ft) 2 

[(2 - ft)\ + 2ft- l]F' + (ft- 1)(F"~ 



2- ft 



V = P', 



H" + 



2-/3 



-H' + XF'H = 0, 



U = V = H = 0, at C = 0, U,V,H -> as C ->■ oo, 



(29) 



(30) 
(31) 
(32) 



which coincides with that in Ref.— except for the fact that we are using the variable 
H = U + W instead of W. This makes sense since the relevant eigenfuncions are such 
that H = O (^e^^- -^ 2 1 2 ^ as £ — > oo (FigfJl right), which is consistent with the asymptotic 

behavior of the variable H (see Figj3] below) as explained above, but not with that of W, 
whose decay is much slower, as e~ g ^^~ a ^ (see eq. (|2"oT) ). 

Now, eq. flBT]) decouples from the remaining three equations and yields the eigenvalue 
A. These eigenvalues will be called Tumin eigenvalues hereafter, and are all positive. The 
two smallest eigenvalues are plotted with solid line and indicated with the subscript T 
in FigfSl left. Note that the first eigenvalue is larger than 1 if < ft < ft c = 0.362, 
which invoking eq. (1271) means that the associated streamwise and wall normal velocity 
components exhibit algebraic growth; the spanwise component instead decays streamwise, 
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as do all velocity components for the remaining eigenvalues. If j3 > /3 C , all eigenvalues are 
larger than one and all flow variables decay streamwise. The eigenfunction components 
U, V, and H associated with the first Tumin eigenvalue are plotted in Figf2j right for 
various representative values of j3. Note that U = —\H\ max U /\U\ ma _ x is always quite close 
to H, which means that U and H are always almost linearly dependent. 

In order to obtain a complete system of eigenf unctions (to, e.g., set all possible initial 
conditions), a second eigenvalue problem must be considered, which results from the 
scaling (cf (ETJ)) 

(U, V, H, P) ~ s 1_A (tf, V, H, p) e -°W«- a '>\ (33) 

Proceeding as above, a second eigenvalue problem results that is a three-dimensional 
extension of its counterpart first considered by Chen and Libby^ in a two-dimensional 
setting. The resulting continuity and spanwise momentum equations must be replaced by 

(1 - X)U + ^ 2 ~_ l f u' + V' = 0, (34) 
, F ~, (g -1)[F+ (2/3 - 1)(F' + Q9 - 1)C 2 F"] - 

+ j^ h+ w^w 

2-/3 2-/3 ' v ; 

but the remaining equations are still eqs. (!29|) and ( |30|) . as are the boundary conditions 
( 132]) . Thus, eqs. (!29|) and ( )34|) are decoupled and provide the eigenvalues, called hereafter 
Chen-Libby eigenvalues. These eigenvalues are all larger than one (the smallest one is 
plotted in Fig. [21 left) and thus they promote streamwise decay; they are also larger than 
the first Tumin eigenvalue, which provides the most dangerous behavior. 

An important property of the two eigenvalue problems considered above is that any 
arbitrary initial condition of the MLBL equations, (Uq, Vq, Hq, Pq), can be expanded in 
terms of the associated modes. In particular, the component of any initial condition on 
the first Tumin eigenmode will be relevant in section IIHI It is given by 



e^P^ Fd ^F'H H T dC 



o 



/ e*=*fi Fdn F'H%.dC 
Jo 



(36) 



as obtained multiplying (with the L 2 inner product) the series expansion by the adjoint of 
the first Tumin mode, Vf, H%, P£) = (0, 0, e ^^ Fdv H T , 0), and taking into account 
that this is orthogonal to the remaining eigemodes. The expression above of the adjoint 
of the first Tumin eigenmode is obtained rewriting eqs. 028p -(j3i p . which provide the first 
Tumin mode, as C-JJt + = 0, C 3 Ht = 0. The adjoint of this problem is 

ciPt = o, c* 2 u* + c;h* = o, (37) 
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where C* stands for the adjoint of the operator Cj. Since A is not an eigenvalue associated 
with the Chen-Libby problem, the first equation of (1371) implies that = 0, which 
means that = _Pf = and is given by C* 3 H^ = 0, where £3 is the adjoint of 
the operator defined by the left hand side of eq. (l3~Tj) . Then, is readily found to be 
H?p = e 4 - 2 < 9 ^° Fc1ii Ht if the usual L2 inner product is used. 

III. STREAMWISE EVOLUTION OF STREAKS FOR < 1 

The numerical method to integrate the MLBL equations (|T3l)-( {TTl) is an extension of 
its counterpart developed in Ref.— . Since the basic steady state converges quite fast to its 
asymptotic value as ( — > 00, the approximation in subsection III Al applies for moderately 
large values of (, say ( > L = 12. In particular, eqs. (l20l and ( 1211) can be used, which 
excluding divergent behaviors as ( — > 00, lead to 

d(V + g(x)V = 0. (38) 

This equation is independent of all MLBL equations. If we substitute the wall normal 
momentum equation by this equation, the resulting system is readily seen to provide the 
correct asymptotic behavior as ( — > 00, analyzed in subsection III Al Thus, we select 
the domain of integration L = 20 > L = 12 and consider eqs. (|T3|) . (|T4|) . and (fl6|) in 
< ( < L, eq. (TL5l) in < ( < L , and eq. (l38l) in L < ( < L. In addition, second order 
spatial derivatives in the resulting system are discretized using second order centered 
differences; first order derivatives of V and P are discretized using second order forward 
differences in eqs.l fTBl) and ( |38i) and second order backward differences in ( fl5|) ; the latter 
are also used to discretize first order derivatives of U and H in all equations. This means 
(noting that both V and P are discretized in Lq < ( < L with forward differences) that 
no boundary conditions for V and P are needed at ( = L\ instead, V(L) and P(L) are 
selected by the numerical code itself. Since both U and H decay extremely fast as ( — > 00, 
the boundary conditions for these at £ = L are U = H = 0. The resulting equations can 
be written as xM(x)<iq/ 'dt = L(x)q, where q = (U,V,H,P) T is the joint flow vector 
and the matrices M and L result from the left and right hand sides of the equations. 
The new logarithmic variable s = log(x/cto) is used to integrate the system from x = Xq 
to x = 1, discretizing s-derivatives with second order forward differences, and using an 
implicit scheme to march in s. 

The resulting numerical scheme is very efficient because the numerical instabilities are 
excluded and the behavior as ( — > 00 is well captured by f )38|) . The scheme is also robust 
because results are insensitive to both the choice of L and L (provided that L > 15, 
Lq > 10, and L is somewhat large compared with L ) and the spatial and s-discretizations 
(provided that the latter be not too fine compared with the former since the mass matrix 
M is singular). Note that having solved the problem for ( — > 00 has allowed to describe 
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accurately the flow very close of the leading edge (x — > 0), without the need of using a 
extremely fine mesh in (, which would have made the method impractical at small x. A 
typical run is carried out in 10~ 9 < x < 1 with Lq = 12 and L = 20, considering 200 
equispaced £-mesh points and performing 400 s-steps, which requires 20 CPU seconds 
using a Fortran90 code in a standard desktop computer. 

A. The most unstable streaky mode for (3 < 1 

Let us define the most unstable streaky mode (MUSM) as that solution of the MLBL 
equations that results from taking as initial conditions at x = Xq <C 1 precisely that 
reconstructed via eq.( T271) using the first Tumin eigenfunction defined by (I28l)- (l32l) . with 
X = Xj. Note that the MUSM is defined up to a constant factor, common to all flow 
variables. 



0.8 
0.6 

\0\ 

0.4 
0.2 
0, 
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FIG. 3. Streamwise, cross flow velocity profiles, and iJ-profile for (3 = 0.353 (that value of (3 con- 
sidered in Ref.—) rescaled with their maxima in < C < oo, at x = 10~ 5 , 10~ 4 , 10" 3 , 10~ 2 , lO" 1 
and 1; arrows indicate increasing values of x. In addition, the asymptotic profiles, reconstructed 
using (|27D and ([28]) -([32]), with A = Xj, for x = 10~ 5 , 10" 4 , 10" 3 , and 10~ 2 are plotted with 
dot-dashed lines using the same rescaling. 

Integration from x = xq = 10 -9 to x — 1 provides the streamwise, cross flow velocities, 
and the if-profile (normalized with their maxima in < £ < oo) illustrated in Fig]3] for 
(3 = 0.353 (that value considered in Ref.—); the results for other values of /3 are completely 
similar. Note that our analysis above captures quite well the right behavior at small x, 
plotted with dot-dashed lines in Figj3] This is possible because of our re-interpretation 
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of Tumin asymptotic results in terms of the new variable H and the use of ([2"T]) . Also 
note that, as occurred in the case of Blasius boundary layer—, the solution exhibits an 
approximate selfsimilarity in the streamwise velocity component and the variable H, which 
(after rescaling with their maxima) are approximately independent of x and approximately 
equal to each other. This illustrates the low- dimensional nature of the MLBL equations 
in the Falkner-Skan boundary layer. 



10 4 
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FIG. 4. Maxima (in < £ < oo) of U, V, W, and H versus x (solid lines, as indicated by the 
labels), for j3 = 0.1,0.2,0.353, and 0.5. The asymptotic behaviors as x — > (see section Bl B|) 
are also plotted with dot-dashed lines. The result of applying random initial conditions, as in 
(|41|) . using the projection on Tumin initial condition in eq. (|36|) . is plotted with dashed lines. 

Figure H] shows with solid lines the maxima (in < £ < oo) of U, V, W, and H vs. x 
for four representative values /3, namely (3 = 0.1, 0.2, 0.353, and 0.5. The above mentioned 
scaling factor (common to all flow variables) has been chosen such that the maximum of 
U vs. x curve in xq < x < 1 is one. The asymptotic behaviors near the leading edge 
W ~ H ~ x^?~ Al and U ~ V ~ x 1 ~ A ^ (see (|27|) ) are also plotted with dot-dashed lines 
for comparison. Note that transient growth is present at /3 = 0.1,0.2, and 0.353, while 
everything decays at /3 = 0.5, as predicted in section III Bl 

The MUSM is somewhat similar to the standard eigenmodes in strictly parallel flows, 
except of course for the final viscous dissipation decay at large x, which is intrinsic to 
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transient growth. In particular, the MUSM is the first of a sequence 

qj = (Uj^H^Pj) for j = l,..., (39) 

which is obtained taking as initial condition at some small xq the Tumin and Chen-Libbby 
eigenmodes associated with the eigenvalues Ai, A2, . . . (sorted together in increasing order). 
Invoking the scalings (J27l) and ( l33l) . such initial conditions are of the form (0, 0, Ht, 0) and 
(Ucl, Vcl, 0, Pcl) for the Tumin and Chen-Libby modes, respectively; in fact, a better 
(but asymptotically equivalent) definition of the initial conditions are obtained using ff271) 
and (l33l) with the remaining components of the Tumin and Chen-Libby modes that have 
been set to zero above. The streamwise evolving modes (139]) satisfy 

IIM > IIQ2II > ■ • ■ inx <x<l, (40) 

provided that Xq *C 1. In addition, any streak can be expanded in series of these modes 
because any initial condition can be expanded in a series of the Tumin and Chen-Libby 
modes. 

All these imply that after a streamwise transient any streak converges to the MUSM, 
namely q x . In addition, the projection of any streak into q : is obtained projecting its 
initial condition into the eigenfunction of the first Tumin mode, as given by eq. fl36l) . This 
is illustrated in FigJH where various streaks are considered that result from imposing at 
various x = x random initial conditions of the type (U, V, H, P) = (Uo, 0, H , 0), with 

Uo = xl^CF" J] ll feC, Ho = x^C'F" £ ll cos k(. (41) 

k=0 k=l 

Here, the coefficients 7^ and 7^ are chosen randomly between 1 and -1. Note also that 
fixing two of the four variables makes sense since the initial condition should satisfy two 
compatibility conditions (one obtained multiplying ffT3"|) by F' and subtracting (TTIj) . and 
the second one, substituting (fT5l) - (|T6l) and the x derivative of (fT4|) into (TT3T) ) . which are 
selected by the equations after a few integration steps if not satisfy initially. These initial 
conditions are applied several times for different values of xq. Results (after rescaling) 
are shown with dashed lines in Fig. HI Note that transients survive until x/xq ~ 10, 
assuming that the initial amplitudes of first and second Tumin modes are comparable. 
This is explained noting that decaying to the MUSM is associated with the second Tumin 
mode (see Figj2j), which compared with the dominant behavior decays as (xq/x) 11 , where 
^ = \t - \t = 0.96, 1.02, 1.13, and 1.25 for = 0.1, 0.2, 0.353, and 0.5, respectively. 

IV. OPTIMAL INTERNAL STREAKS FOR (5 < 1 

Optimal streaks are those streaks that show the maximum amplification. Finding 
optimal streaks is made maximizing the perturbed kinetic energy gain G = E m /E out 
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between a generic section x 



X ; 



and x — 1, which invoking §5§ is written as 



G,, 



max 





J 


poo 

I U 2 dy 




x=l 







z?n J 


poo 

I (y 2 + W 2 )dy 







(42) 



A factor Re is omitted in the numerator, where a 0(1/Re) term depending on the cross 
flow velocities is also neglected. The latter cannot be done in the denominator because, 
in fact, the 0(i?e)-term becomes negligible in the maximizers. Thus, following Luchini^, 
Tumin and Ashphis^ set to zero that term imposing the additional condition that U = 
at x in . Here instead, we shall retain the 0(-Re)-term in the denominator, checking (as in 
Ref.— ) that results are independent of Re provided that Re be large. 

Note that maximizers of (1421) depend on x in , with the asymptotic result as x- lTL — > 
being most relevant. Invoking the scaling (fTTI) and the definition fl26l) . eq. fT42|) is rewritten 
as 

U 2 d( 



x 



g+i 

2-/3 



G, 



la 



max-F 



Uo 



2-/3 



X 



2 

2-/3 



Re 



U 2 d(+ / [V 2 + (H - V) 2 }d( 



(43) 



2-/3 
x=a 1 ~P c 



'0 JO 

This quotient could be maximized with the method used in Refs.— _ — , who wrote the 
quotient as (qout'Qout'Qout) / (Qin'Qm'Qm)' where q is the flow state vector q = (U, V, H, P). 
Using the action U associated with the streamwise evolution dynamical system associated 
with the MLBL equations, defined as q out = U-q in , and its upstreamwise-evolving adjoint, 
defined such that (U ■ q in ) T = qj n ■ U*, the quotient ( l4~2l) is rewritten as G = (qj n ■ U* • 
Qout ' U ■ q in ) / (q[ n ■ Q in ■ q in ). Maximizing this leads to the generalized eigenvalue problem 
U* • Q out ■ U ■ q in = GQ in • q in , whose maximum eigenvalue provides the maximum gain 
and the associated eigenf unctions, the maximizers. This problem is iteratively solved as 
q inn+1 = • U* ■ Q out ■ U ■ q inn , which converges quite fast provided that the two first 
eigenvalues of the above mentioned generalized eigenvalue problem are not too close each 
other. 



A. A simple optimization method to calculate optimal streaks 

The main difficulty in the method outlined above is the need of using the adjoint of 
the MLBL equations, which was made in Refs.— _ — for the original linearized boundary 
layer equations, written in terms of the unsealed wall normal coordinate y. A similar 
adjoint equations could be derived for the MLBL equations in this paper, but doing that 
in an efficient way would require to analyze the asymptotic behaviors of these equations 
as C — > oo and 0. Instead, we propose here a simpler method that is new in 

this context to our knowledge. This method is based on the observation above that the 
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infinitely many solutions of the MLBL equations can be classified as a series of streamwise 
evolving 'modes' indicated in ( 139]) and that they satisfy ( l4"U|) . Thus, if a particular solution 
is written as linear combination of the first n such solutions, the resulting error scales with 
(x; n /x) An//Al , which becomes exponentially small as x/x- m is moderately large if A n /Ai is 
moderately large. This means that higher order modes can be safely neglected when 
optimizing (143j) because these will contribute to increase the denominator of (1431) but will 
have a negligible effect in the numerator. Now, the sequence Ai, A 2 , . . . increases rapidly, 
meaning that retaining a few modes would be enough to define well generic solutions of 
the MLBL equations. Thus, the method we propose to maximize (143]) is as follows: 

Step 1. Take the first n Tumin and Chen-Libby modes, ordered such that the associated 
eigenvalues are sorted in increasing order, as explained right after eq.( |39l) . 

Step 2. Calculate the n solutions of the MLBL equations, 

(U h Vl, H h Pi),..., {tin, V n , H n , P n ). (44) 

2-/3 

obtained taking as initial conditions at x — x- m = a 1 -> 3 Xi n the n modes calculated 
in step 1. 

Step 3. Replace the expansion 

n 

(U, V, H,P) = J2 a i(U S , % Hj, Pj), (45) 

3=1 

into ( 1431) . to rewrite the quotient appearing in this equation as 

S+l 

G max («) = Y k 2 in Jk J \ , (46) 

E2, =1 (x£'1teE$ + E$ )a jak 



where E°f = UjU k dC at x = , Ef k x = J °° UjU k d( at x = a~^x- m , and 
E fk = iriWk + (Hj - V 3 )(H k - V k )]d( at x = a^x in . 

Step 4. Equation f)46p is the ratio of two n-th order quadratic forms, which is maximized as 
usually, solving the generalized eigenvalue problem 



i+l 
2-13 



n n , 2 «. 

J^E^W = GJ2 [<n^Ef + Ef\a k , (47) 

k=l k=l ^ 



whose maximum eigenvalue G max provides the maximum gain in (143]) ; the eigenvec- 
tor yields the maximizer of (1431) using ( 145]) . 
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FIG. 5. Rescaled (with the maximum in a) maximun perturbed energy gain (|43j) for (3 = 0.1 
(top, left), 0.2 (top, right), 0.353 (bottom, left), and 0.5 (bottom, right), taking the initial stage 
at Xi n = 10~ 5 (dot-dashed lines), 10 -2 (dashed lines), and 0.25 (solid lines). Results from Ref^ 
for (3 = 0.353 and x m = 0.25 are plotted with plain circles. 

The method is tested considering the cases j3 = 0.1,0.2,0.353 (the case considered in 
Ref.—), and 0.5, with three values of the initial streamwise stage, namely x in = 10~ 5 , 10~ 2 , 
and 0.25 (as in Ref.—). The associated rescaled maximum perturbed energy gain curves 
obtained applying the method described above are given in FigfSJ in order to facilitate 
comparison for the various values of x- m , the energy gain is rescaled with its optimal value 
(plotted vs. (3 in Figj6]below). One, three, and five modes are enough for x- m = 10 -5 , 10 -2 , 
and 0.25, respectively, at the four considered values of /3; retaining more modes provides 
results that are plot indistinguishable. As a reference, the counterpart obtained in Ref.— 
at j3 = 0.353 and x in = 0.25 are plotted with plain circles; in fact, Tumin and Ashpis took 
x; n = 0.111 and x out = 0.444, which must be compared with x- m = 0.25 and x out = 1 with 
the adimensionalization in the present paper. 

This method provides the maximum gain in a quite fast and robust way. Note that the 
method does not require any calculation of adjoint equations. Instead, only a few Tumin 
and Chen-Libby eigenmodes are needed. But these are not really necessary, noting that 
the n solutions of the MLBL equations considered in step 2 can be replaced in the method 
by any set of n solutions with initial conditions at x = x- m = a^x^ that are linearly 
independent and exhibit a significant projection into the n first Tumin and Chen-Libby 
eigenmodes. For instance, we can integrate the MLBL equations taking at some smaller 
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value of x (say, /10) the following n initial conditions (cf (14"T|) ) 



(U k ,V k ,H k ,P k ) 



x 



1--^ta2 t?h ^„ us n ,5,2-/3 A T / .2 I? // 



X F sin fcC,0,Xo 



( F cosfcCO 



(48) 



for k = l,...,n, and considering in step 2 these solutions for x > x- m . The resulting 
modified method, retaining the same numbers of modes as above produce the same results 
in Figj5] to plot accuracy. Since this modification provides a simpler method, it is the 
resulting modified method that we propose in this paper. 

Note that the number of required modes in Figj5] decreases as x- in decreases. In par- 
ticular, just one mode is enough if x- m < 10~ 5 , and furthermore, the resulting maximizer 
coincides with the MUSM defined in section UlI Al as was to be expected. In fact, taking 
into account the behavior of the MUSM as x — > (given by the eigenfunction associated 



with the first Tumin eigenvalue), the following asymptotic value of the gain 
obtained that does not require any maximization process 



la 



2-/3 


" POO 

/ u 2 dC 

Jo 


2-/3 
x=a 1 ~P 




" POD 
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H 2 d( 




2-3 











X:. 



<3-l 
"2-/3 



as X\ 



0. 



is readily 



(49) 



2-fl 

Here, U and H are the associated components of the MUSM, at x = a 1 -/ 3 and x = a 1 ~i 3 Xi n , 
respectively. The indicated asymptotic behavior results from the scaling (I2j 

of the first Tumin mode, associated with the eigenvalue A^. 



2-/3 
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FIG. 6. Optimal spanwise wavenumber a opt vs. /3 at x; n = 10 5 (dot-dashed line), 10 2 (dashed 
line), and 0.25 (solid line). Plain circles at (3 = 1 result from calculations in the Appendix. 

To complete the results above, the optimal gain and the associated value of the spanwise 
wavenumber are plotted vs. /3 in Figj6]for x- m = 10~ 3 , 10~ 2 , and 0.25. The limiting values 
at (3 = 1 result from the analysis in the Appendix, at the end of the paper. Note that 
a opt depends non- monotonously both on (3 and on x\ n . 

Summarizing, a simple optimization method has been proposed to calculate optimal 
streaks that does not rely on the adjoint formulation. Instead, minimization is made 
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on the low dimensional manifold spanned by a few solutions of the direct problem. In 
fact, five solutions are enough in the streaks of the Falkner-Skan boundary layer for 

l£ 0.25. In addition, only one such solution is enough for small x- m (x- m 
meaning that no optimization process is necessary in this limit. We believe that the 
method proposed above is also useful to treat related transient growth problems in Fluid 
Dynamics, since the main ingredient that allowed constructing such method is usual in 
these problems. This ingredient is that the behavior in the growth stage of transient 
growth are described by an eigenvalue problem (Tumin problem in the case considered 
in this paper), whose eigenvalues are somewhat separated. And furthermore, the method 
could work in nonlinear, time dependent parabolic problems as well, since the dynamics 
of parabolic problems is low dimensional at large time^. 

V. CONCLUDING REMARKS 

An analysis of optimal streaks in the Falkner-Skan boundary layer has been performed 
concentrating in various issues that are now summarized, and are expected to apply to 
related transient growth boundary layer problems: 

• The careful analysis of the behavior of streaks near the free stream was necessary 
to construct a quite efficient numerical scheme, which allowed for integrating from 
extremely small values of the streamwise coordinate. The analysis of the free stream 
behavior also allowed for making the correct interpretation of the behavior near the 
leading edge, in terms of a new variable H that behaves as the streamwise vorticity 
near the leading edge. 

• The already known approximate selfsimilarity of the solution has been completed 
in terms of the new variable H. 

• Streaks behave as 'modes', which can be classified according to their behavior near 
the leading edge. The analogy with standard modes in strictly parallel flows becomes 
clear comparing the cases j3 < 1 (considered in most part of the paper) and (3 = 1 
(considered in the Appendix). The most dangerous mode was called the MUSM 
and played an essential role in understanding the streamwise evolution of streaks. 

• A simple optimization method has been proposed that does not rely on the adjoint 
equations and allows for the fast and precise computation of optimal streaks. 

• Optimal streaks can be directly defined from the MUSM, without the need of any 
optimization process, if the initial stage is sufficiently close to the leading edge. 

• All the above is somehow a pre-requisite to derive simple, yet sufficiently precise de- 
scriptions of the interaction between longitudinal streaks and transversal Tollmien- 
Schlichting modes, which can be still considered a major open problem in the field. 
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• The analysis above applies to internal streaks, resulting from, e.g., obstacles near 
the leading edge. External streaks forced by perturbations in the outer stream show 
a different behavior near the leading edge and require a different treatment, which 
is currently under research. 
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APPENDIX: THE CASE > 1 

In the case (3 = 1 the spanwise wavenumber a cannot be eliminated from the LBL 
equations (|6])- ([9]), but the right hand sides of these are independent of x and thus can be 
expanded in terms of normal modes as (cf eq. ()27j) ) 



(1993). 



oo 




(50) 
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where c.c. stands for the complex conjugate and the modes and exponents are given by 
the following eigenvalue problem 

(1- X)U + V'-aW = 0, (51) 

U" + FU' + [(A - 2)F' - a 2 ]U - F"V = 0, (52) 

V" + FV' + [(A + 1)F' - a 2 }V = P', (53) 

W" + FW 1 + (XF' - a 2 )W = aP, (54) 

U = V = W = at C = and oo, P = at ( = oo. (55) 

Note that as a — > two scalings are possible. Rescaling W as W/a and setting a = 0, 
(I5TT) and (1541) lead to the Tumin eigenvalue problem (|28|) - (|3T|) . while setting a = in 
(ED) and ([55]), the Chen-Libby eigenvalue problem pfll-CT, and (155|) results. This 
means that the limiting values of the eigenvalues of ()5ip and fl55l) as a -)• are Tumin 
and Chen-Libby eigenvalues. 




FIG. 7. The real (left) and imaginary (right) parts of the first three eigenvalues of (|51[) - fj55() vs. 
a. Plain circles at a = correspond to the Tumin and Chen-Libby eigenvalues at (3 = 1, plotted 
with plain circles in Figf2J 

The eigenvalues of (IBTj) and (1551) are generally complex, as shown in FigJTJ where the 
real and imaginary parts of first three eigenvalues are plotted vs. a; note that the second 
and third eigenvalues are complex in the intervals 0.16 < a < 0.96 and 0.27 < a < 0.94, 
respectively. 

The limiting values as /3 — > 1 of the optimal wavenumbers in Figj6] are calculated 
minimizing the counterpart of the ratio (f46j) . using n of the modes ([5"0]) . which is 

^ I \ — Zlj,fc=l x inEjk t °'j a k , s 

where overbar stands for the complex conjugate, = UjUkdC, at x = a, E 1 ^ 1 = 
f °° UjUkdC at x = ax in , and E^ 2 = / °°[t^Vfc + (Pj - Vj)(H k - V k )]d( at x = ax in . The 

21 



Optimal streaks 

resulting optimal gain calculated from the generalized eigenvalue problem (cf (117)1 ) 

n n 

E = G E + Efk 2 ) o*, (57) 

fe=i fc=i 

whose maximum eigenvalue G max provides the maximum of the energy gain ( 1571) . Note 
that although the amplitudes are complex, the eigenvalues are real. This is because 
the matrices appearing in (1571) are Hermitian. Using these, the maximum perturbed 
energy gain is calculated from the initial stages x- m = 10~ 5 , 10~ 2 , and 0.25, which peak 
at a = 0.38, 0.45, and 0.60, respectively. These are precisely the limiting values of their 
counterparts calculated for f3 < 1 in subsection IIV Al and plotted in FigJHJ 

Even though the case j3 > 1 does not produce transient growth, a brief description is 
given here to complete the analysis. This case yields a different behavior of the solutions 
of the MLBL equations (fT3]) - (fT7|) as x — > 0. This is because g(x) = x 1 ^ — » oo as x — > in 
this case. Since g(x) is large, U, V, H, and P converge to zero exponentially fast as g(x)( 
is large, which means that activity is concentrated in a thin wall-normal layer near the 
wall. The x-behavior in this stage is given by an eigenvalue problem that is omitted here 
for the sake of brevity. As ct — > oo instead g(x) is small and the Tumin and Chen-Libby 
eigenvalue problems are recovered. Thus, the solution converges to the first Tumin mode 
as x — > oo. In order to illustrate all these, the MLBL equations are integrated for (3 = 1.5, 
between x = Xq = 0.1 and 10 3 . The maxima of \U\, \V\, \H\, and \W\ are plotted with 
solid lines in FigJHlleft. Note that all these decay quite fast in an initial stage, followed by 
a transition to the behavior given by the first Tumin eigenmode, plotted with dot-dashed 
lines. The irregular behavior at small x is not a numerical artifact, but is due to the fact 
that we are plotting maximum values, and the maxima alternated between two positions, 
as illustrated in the right plot, where the rescaled profiles of \U\ are plotted for various 
values of X, clS indicated. 
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FIG. 8. The case (3 = 1.5. Left: Maxima (in < ( < oo of \U\, \V\, \H\, and \W\ vs. x. 
The upper and lower straight dot-dashed lines provide the asymptotic behaviors of \H\ and \U\, 
which invoking (llip . fjl2|> . and (|27|) are and x 1 ~ x /g(x), respectively. Right: A sequence of 
profiles of \U\ at nine equispaced values of x between x = 0.1389 and 0.1995, x\, . . . ,x$. Due 
to the exponential decay, these profiles cannot be plotted in the same plot. Instead, two of 
them are given in each plot, which also contains the last profile in the former subplot, namely 
the four subplots provide the profiles at x = xi,X2,x% (top, left), x = £3, £4, £5 (top, right), 
x = x§,xq,x 7 (bottom, left), and x = xj,x 8 ,xg (top, right). Arrows indicate increasing values 
of x. 
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